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Abstract 

Self-seeding schemes, consisting of two undulators with a monochromator in be- 
tween, aim at reducing the bandwidth of SASE X-ray FELs. We recently proposed 
to use a new method of monochromatization exploiting a single crystal in Bragg- 
transmission geometry for self-seeding in the hard X-ray range. Here we consider 
a possible extension of this method to the soft X-ray range using a cell filled with 
resonantly absorbing gas as monochromator. The transmittance spectrum in the 
gas exhibits an absorbing resonance with narrow bandwidth. Then, similarly to the 
hard X-ray case, the temporal waveform of the transmitted radiation pulse is char- 
acterized by a long monochromatic wake. In fact, the PEL pulse forces the gas atoms 
to oscillate in a way consistent with a forward-propagating, monochromatic radia- 
tion beam. The radiation power within this wake is much larger than the equivalent 
shot noise power in the electron bimch. Further on, the monochromatic wake of 
the radiation pulse is combined with the delayed electron bunch and amplified in 
the second undulator. The proposed setup is extremely simple, and composed of as 
few as two simple elements. These are the gas cell, to be filled with noble gas, and 
a short magnetic chicane. The installation of the magnetic chicane does not perturb 
the undulator focusing system and does not interfere with the baseline mode of 
operation. In this paper we assess the features of gas monochromator based on the 
use of He and Ne. We analyze the processes in the monochromator gas cell and out- 
side it, touching upon the performance of the differential pumping system as well. 
We study the feasibility of using the proposed self-seeding technique to generate 
narrow bandwidth soft X-ray radiation in the LCLS-II soft X-ray beam line. We 
present conceptual design, technical implementation and expected performances 
of the gas monochromator self-seeding scheme. 
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1 Introduction 



Self-seeding schemes have been studied to reduce the bandwidth of SASE X- 
ray PEL lUt] - 10] • In general, a self-seeded EEL consists of two undulators with 
an monochromator located between them. The first undulator operates in 
the high gain linear regime starting from the shot noise in the electron beam. 
After the first undulator, the output radiation passes through the monochro- 
mator, which reduces the bandwidth to the desired value, smaller than the 
EEL bandwidth. While the radiation is sent through the monochromator, 
the electron beam passes through a bypass, which removes the electron mi- 
crobunching introduced in the first undulator and compensates for the path 
delay created during the passage in the monochromator. At the entrance 
of the second undulator, the monochromatic radiation pulse is recombined 
with the electron beam, and amplified up to saturation. The radiation power 
at the entrance of the second undulator is dominant over the equivalent shot 
noise power, so that the bandwidth of the input signal is smaller than the 
bandwidth of EEL amplifier. 

The self-seeding scheme proposed for the VUV and soft X-ray region makes 
use of a grating monochromator |jlt] • The presence of such monochromator 
introduces a path delay with respect to the straight path, which has to be 
compensated with the introduction of a bypass of the length of about 30 
m for the case of the LCLS-II [5]. This requires relevant modifications of 
the baseline. We recently proposed to use a new method of monochrom- 
atization exploiting a single crystal in Bragg-transmission geometry for 
self-seeding in the hard X-ray range ||6(l-[|10l]. A great advantage of such 
transmission geometry method is that it introduces no path-delay of X-rays 
in the monochromator, thus avoiding the need for a long electron bypass. 

All VUV and soft X-ray monochromators operate in the frequency domain 
as bandpass filters. In this paper we propose, instead, a new method of 
monochromatization based on the use a cell containing resonantly absorb- 
ing gas in the transmission direction. Such cell works, in analogy with 
the single-crystal monochromator method for the hard X-ray, as a band- 
stop filter. Then, such problems as extra path-delay of radiation pulse, heat 
loading and monochromator alignment are solved. The proposed setup is 
extremely simple, and composed of two simple elements: a gas cell and 
a short magnetic chicane, Eig. 1. As for the hard X-ray case, the magnetic 
chicane accomplishes three tasks. It creates an offset for gas cell installation, 
it removes the electron micro bunching produced in the first undulator, and 
it acts as a delay line for the electron bunch. Thus, using a cell with rare gas 
installed within a short magnetic chicane in the baseline soft X-ray undula- 
tor as a bandstop filter, it is possible to decrease the bandwidth of the VUV 
radiation down to the Eourier limit. 
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Fig. 1. Design of the VUV - soft X-ray FEL undulator system for the narrow band- 
width mode of operation. The scheme is based on the use of the self-seeding scheme 
with a gas monochromator. The magnetic chicane absolves three tasks. It creates an 
offset for installing the gas monochromator, it removes the electron micro-bunching 
produced in the first undulator, and it shifts the electron bunch onto the monochro- 
matic wake. After the gas cell, the monochromatic wake of the radiation pulse is 
combined with the delayed electron bunch, and amplified in the downstream un- 
dulator. The downstream undulator itself consists of two parts: a first part, in which 
the monochromatic seed signal is amplified in the linear regime and a second part, 
in which the harmonic of the seed signal produced in the first part is amplified up 
to saturation. 

We will show that the proposed technique can be implemented in the soft 
X-ray region, exploiting the tunable-gap soft X-ray baseline undulator. We 
propose to first perform monochromatization at 20 - 40 nm with the help 
of the gas-cell monochromator, and subsequently amplify the radiation in 
the first part the of output undulator. The amplification process can be 
stopped at some position well before the FEL reaches saturation, where the 
electron beam gets considerable bunching at the 2nd, 3rd and 4th harmonic 
of the coherent radiation. The undulator downstream that position can then 
be tuned to the harmonic frequency to amplify the radiation further to 
saturation. Fig. 1. 

In this paper we present a detailed study for applying the self-seeding 
scheme with gas monochromator to the LCLS-II setup, in order to generate 
narrow bandwidth, soft X-ray radiation in the proposed soft X-ray beam 
line. The tentative design for the proposed technique at LCLS-II allows to 
generate fully coherent radiation at a wavelength of about 5 nm. We focus on 
the conceptual design, on the technical implementation and on the expected 
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Fig. 2. Gas cell with He, Ne and Ar as working gas, working as a bandstop filter 
for the transmitted VUV SASE radiation pulse. 

performances of the self -seeding scheme using the gas monochromator. 



2 Principles of the self-seeding technique based on the use of a gas 
monochromator 



In this section we illustrate more in detail our novel method of monochrom- 
atization, based on the use of a cell containing resonantly absorbing gas. As 
already said, this technique takes advantage of the transmission geometry, 
where no extra path-delay for the radiation pulse is present. The principle 
of the new method of monochromatization is very simple and is illustrated 
in Fig. 2 and Fig. 3. An incident SASE pulse coming from the first undulator 
impinges on the gas cell. When the spectral contents of the incoming SASE 
beam satisfies the photo-absorption resonance condition of the gas, the gas 
cell actually operates as a bandstop filter for the transmitted VUV SASE 
radiation pulse, see Fig. 3. 

Obviously, if we use a bandstop filter there is no monochromatization in 
the frequency domain. However, the temporal waveform of the transmitted 
radiation pulse shows a long monochromatic wake. The duration of this 
wake is inversely proportional to the bandwidth of the absorption line in the 
transmittance spectrum. It is then possible to reach a bandwidth limited seed 
signal by exploiting a temporal windowing technique, as indicated in Fig. 
3. Since we deal with a high-gain parametric amplifier where the properties 
of the active medium, i.e. electron beam, depend on time, the temporal 
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Fig. 3. Forward radiation in the gas cell. The transmittance spectrum in the gas cell 
shows an absorption resonance with narrow (4-10 meV) line width. When the spec- 
tral contents of the incoming radiation beam overlap the absorption resonance, the 
temporal waveform of the transmitted radiation pulse exhibits a long monochro- 
matic wake. The duration of the wake is inversely proportional to the bandwidth 
of the absorption resonance. 

windowing concept can be practically implemented in a simple way by 
delaying the electron bunch at the position where the frequency spectrum 
of the transmitted pulse experiences a strong temporal separation. In other 
words, the magnetic chicane in Fig. 1 shifts the electron bunch on top of the 
monochromatic wake created by the bandstop filter. By this, it is possible to 
seed the electron bunch with radiation pulse characterized by bandwidth 
much narrower than the natural FEL bandwidth. The concept is similar to 
that developed in |@]-|10 ] for the hard X-ray case. However, in the proposed 
scheme no optical elements are used, and problems with alignment or heat 
loading do not exist at all. While the radiation is sent through the gas cell, 
the electron beam passes through the magnetic chicane, which accomplishes 
three tasks by itself: it creates an offset for the gas cell installation, it removes 
the electron micro bunching produced in the first (SASE) undulator and, as 
already said, it acts as a delay line for the implementation of the temporal 
windowing. 

For the VUV wavelength range, a momentum compaction factor J^sg in the 
order of 100 y.vs\ is sufficient to remove the microbunching in the electron 
beam. As a result, the choice of the strength of the magnetic chicane only 
depends on the delay that we need to introduce between electrons and 
radiation. In our case, this amounts to 100 /.im, see Fig. 1. Such delay is 
small enough to be generated by a short 3 - 4 m long chicane to be installed 
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in place of a single undulator module at a facility like the LCLS-II. Such 
chicane is, however, strong enough to create a sufficiently large transverse 
offset of 1 - 2 cm for installing the gas cell. 



3 VUV Photo-absorption spectra of the rare gases 



In this article we will take advantage of autoionizing resonances in rare 
gases. We consider He, Ne and Ar as possible choices, presenting cross 
sections and tables of line parameters. Additionally we will discuss about 
Doppler broadening as an effect, which can be accounted for in this most 
general level of descriptioriH 

The phenomenon of autoionization is well-known in literature (see for ex- 



ample textbooks like IHO, |12|1 and references therein, just as examples of 
reviews). Let us consider the Helium atom as a case study. The interaction 
with a soft X-ray electromagnetic pulse has the effect of ionizing the atom. A 
first ionization threshold, where the remaining electron in the ion is found in 
the state with principal quantum number N = 1, is found at 24.6 eV, while 
a second ionization threshold with N = 2 is found at 65.4 eV. Therefore, 
after interaction with a photon with energy between these two thresholds, 
one finds the Helium ion in the ground state Is, and a free electron with 
kinetic energy equal to the photon energy diminished of 24.6 eV. This single 
open channel can be reached either directly or via decay of a doubly excited 
state, called autoionizing state. A set of autoionizing states can be excited at 
different energies. If the photon energy increase between 65.4 eV to 72.9 eV, 
not only the previously discussed open channel is possible. In fact, above 
the second ionization threshold a second open channel is present, charac- 
terized by the Helium ion with the bound electron with N = 2 (2s state), and 
free electron energy equal to the photon energy diminished by 65.4 eV. In 
this photon energy range we have a set of autoionizing resonances,which 
interfere with two open channels (Is or 2s). More in general, as the pho- 
ton energy increases below the double ionization threshold at 79 eV, other 
channels corresponding to the Helium ion with higher numbers of N (Nth 
ionization threshold) can be reached. 

The autoionizing He* state corresponds to a double electronic excitation 
where the two electron principal numbers become Ui = N and n2 = N,N + 
1, with N > 1. Each choice of N = 2,3... corresponds to a set of Rydberg 
series. Dipole transitions allow 2N + 1 series of ^Pq autoionizing He* states. 



^ The Doppler broadening depends on the gas temperature only (room tempera- 
ture in our case) and not on density or pressure, which will be specified only later, 
when discussing the differential pumping system design. 
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These series of resonances admit as a limiting energy the Nth ionization 
threshold corresponding to the an excited state with quantum number N 
for the He'^ ion. Since for a given choice N series admits the Nth ionization 
threshold as a limit, this means that any autoionizing He* state decays to the 
He* state with quantum number N - 1. 

In this paper we will focus on the series of autoionizing resonances converg- 
ing to the N = 2 levels of Helium ion. This means that we have 2N + 1 = 3 
Rydberg series, all converging to the second ionization threshold, corre- 
sponding to 65.4 eV. These three series all begin with the two Helium 
electrons characterized by ni =2 and n2 = 2, which corresponds to an 
energy of 60 eV. There are three possible choices for the {li,l2) quantum 
numbers: (2s, np) ^Pq, (2p, ns) ^Pq, and (2p, nd) ^Pq. However, the levels of the 
(2s, np) ^Pq and (2p, ns) ^Pq wave functions are nearly degenerate, and the 
doubly-excited states are better described in terms of constructive and de- 
structive superpositions of these wave functions, finally yielding the three 
Rydberg series: (sp,2n+) ^Pq, {sp,2n-) ^Pq, and {2p,nd) ^Pq. Since these three 
series converge (but no reach) the second ionization threshold, all the cor- 
respondent doubly excited states He* decay to an Helium ion in the ground 
state. Out of these three series, only the first (sp,2n+) will be of interest to 
us, as the others are too weak. 



3.1 Photo-absorption cross-section of He for photon energies 60 eV to 65 eV 



As we just said, autoionizing resonances results from the decay of doubly 
excited Rydberg states of He* into the continuum, i.e. He* — > He*e~ . Since 
the continuum can also be reached by direct photoionization, both paths add 
coherently, giving rise to an interference. This interference is related with 
the typical Fano line shape for the cross-section as a function of energy [13]. 
A more detailed physical explanation of the entire process of interaction of 
the PEL beam with the gas is given in Section 6 and 7. Here we limit to report 
the cross section for these series, which can be modeled by the expression 



o{K,q,Y) = cTf,(A) 



(1) 



where the energy-dependentbackground cross-section expressed in Megabarn 
{1Mb = 10"^^cm^) is given by 



(7fc(A) = -0.05504 - 1.3624 ■ lO'^A + 3.3822 • 10"^ , (2) 
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Table 1 

Parameters determining the Fano autoionization profiles for the {sp, 2n+) autoion- 
izing series of Helium, following [14]. 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


r(eV) 


0.0378 


0.0083 


0.0038 


0.0014 


0.0008 


0.0005 


0.0003 


0.0002 


0.0001 


Eres (eV) 


60.14 


63.655 


64.466 


64.816 


64.999 


65.108 


65.181 


65.229 


65.263 


q 


-2.57 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 


-2.5 
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Fig. 4. Fano profiles for the (sp. 2n+) autoionizing series of Helium. The cross-sec- 
tions are calculated following |l4ll . 



with A the radiation wavelength in Angstrom units, while the reduced 
energy &„ is defined as 



£„ — 



2{ER„-hc/A) 



(3) 



The asymmetry index qn, the energy of the nth resonance Er„, the resonance 
width r,„ have been the subject of several calculations. 

The cross-section for the series {sp,2n+) can be calculated using Eq. (1), 
with the parameters given in Table 1, reproduced from [14]. Also, the pho- 
toabsorption cross section a is linked to the light attenuation through a gas 
medium of column density HqI, I being the length of the cell and Uq the gas 
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density. This allows cross-section measurements via the Beer-Lambert law. 



I{a)) = Ioexp[-nolo{a))l (4) 

where Iq is the incident intensity, and I{oj) is the attenuated intensity of the 
transmitted light at frequency co. The values of parameters obtained from 



experiments are in general consistent with calculations these uM- 



The resulting cross-section for the {sp,2n+) autoionizing series of Helium 
is shown in Fig. 4 as a function of the energy. The availability of a series 
resonances allow for tunability of the scheme. However, not all resonances 
can be used. In fact, they should be located further away from each other 
than the PEL bandwidth, and be narrow enough to guarantee the application 
of the temporal windowing process. In our case, a rule of thumb requires 
a width narrower than 10 meV. In this case, only the second and the third 
Helium lines are suitable, as the first one is too wide, and starting from the 
fourth there is not enough separatiorlZl 



3.2 Photo-absorption cross-section ofNefor photon energies 45-49 eV 



In this paper we also study of possibility to use the autoionizing resonances 
of Neon from 2s^2p^ to Islp^np , with n = 3 - 5 This allows our setup 
for further tunability. It should be noted that there is a difference in the 
calculation of the cross-section, when compared with the Helium case for 
energies below the second ionization threshold. In fact, in that case only 
one channel is open, meaning that the final state is the Helium ion in the 
ground state. The two paths towards this channels interfere and one is left 
with the Fano line considered before. However, when the energy increases, 
or one considers another gas like Neon, other channels open up. Then, as 



explained for example in |16|1 , transitions to the continuum may or may not 
interact with the discrete autoionization state, and the expression for the 
cross-section has to be modified as: 



O = Ob ^ ^ + Oa , (5) 

where Ob and Oa are the cross-sections of the interacting and non-interacting 
transitions to the continuum. Defining the total cross-section oq = Oa + Oh, 



^ Such smaller separation maybe useful at VUV facilities, in order to seed at mul- 
tiple photon wavelengths with separation of a fraction of eV, similarly to what has 



been proposed in ilOfl. 
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Table 2 

Parameters determining the Fano autoionization profiles for the {sp, 2n+) autoion- 
izing series of Helium, following [14]. 



1 L 




4 




oo (Mb) 


8.6 


8.0 


8.2 


r(eV) 


0.013 


0.0045 


0.0002 


Eres (eV) 


45.546 


47.121 


47.692 


q 


-1.6 


-1.6 


-1.6 




0.7 


0.7 


0.7 



which one would measure without autoionization, Eq. (5) can be cast in the 
form 



where we also defined the correlation coefficient (equal to unity for the 
Helium resonances under study): 



The relevant resonances can be calculated using Eq. (27), with the param- 
eters given in Table 2, reproduced from [jisl] . All three resonances look 
interesting for the implementation of our technique, although the one cor- 
responding to n = 3 is somewhat larger than 10 meV. It should also be 
remarked that the presence of a correlation parameter different from unity 
leads to a fundamental difference in the Fano profiles, that is, the cross- 
section never goes to zero, at variance with the Helium case. 



3.3 Photo-absorption cross-section ofArfor photon energy equal to 28.5 eV 



Another possible choice |16|1 , suitable for our purposes would be the res- 
onance of the Argon gas at 28.5 eV. Such resonance can also be calculated 
with the help of Eq. (27), and is characterized by = 0.85, q = 0.17 and has 
a width r = 12.6 meV. 
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3.4 Doppler broadening 



To conclude the present Section we discuss Doppler broadening of the Fano 
lines. Due to the finite temperature T of the gas, the frequency distribution 
of an ensemble of radiating atoms will experience a spread [17]. Assuming 
the gas ensemble in thermal equilibrium, the Doppler shift is linked to the 
Maxwell-Boltzman probability distribution of velocities of the atoms. The 
distribution of frequencies turns out to be a Gaussian, with a peak at the 
resonance frequency and a FWHM given by 



where Eq is photon energy at resonance, T is the temperature measured in 
Kelvin degrees, and M is the molecular weight. Assuming a temperature 
of 300 K, the Doppler broadening is pressure independent and turns out to 
amount to about 0.4 meV for He ( molecular weight M = 4 ). For resonances 
under study, Fois substantially smaller than the natural width F ~ 4 - 10 
meV. It follows that simulations directly performed from the calculated Fano 
profiles constitutes a good approximation of the exact profile, which should 
be obtained by convolving the Fano natural line with the Gaussian profile 
due to Doppler shift. 



4 Transmittance of the gas cell used in monochromatization 

As discussed in Section 3, if a monochromatic electromagnetic pulse of 
intensity Jq arid frequency co impinges on a cell of length I, filled with a gas 
with density no, the transmitted intensity obeys the Beer-Lambert law, Eq. 
4. As a result, by comparison with Eq. 4, the modulus of the transmissivity 
can be defined as 



where the shape of o{co), which follows a Fano profile, has been discussed 
in the previous Section, and the product n^l is known as the column-density 
of the gas cell. In the following we will choose a column-density equal to 
nol = lO^^cm"^, and we will restrict ourselves to the third (n = 4) resonant 
line of the {sp, 2n+) ^Pq Rydberg series for Helium, whose cross-sections can 
be calculated with the help of Eq. (1) and Eq. (2). 

Once |T| is known, it is possible to use Kramers-Kroning relations to recover 
the phase. In fact one can write 



Tu = 7■W-'Eo^Jf/M, 



(8) 



\T\ = exp[-no/cT/2] , 



(9) 



12 




19.220 19.225 19.230 19.235 19.240 19.245 

A.[nm] 



Fig. 5. Modulus and phase of the transmissivity of Helium around the n = 4 line 
of the {sp,2n+) ^Pq Rydberg series. The modulus has been calculated according to 
Eq. (9), while the phase is recovered with the help of the Kramers-Kroning relation 
according to Eq. (11). 

ln[r(a;)] = ln[|T(a;)|] + iO{a)) = -nlo/2 + zO(a;) . (10) 



Note that T{a)) = T{-a)) implies that \T{a))\ = \T{-a))\ and that O{co) = 
-0{-oj). Therefore, using Eq. (10) one also has that ln[T(a;)]* = ln[r(-a;)]. 
Then, application of Titchmarsh theorem shows that the analyticity of 
ln[|T(Q)|] on the upper complex Q-plane implies that 



n J co'^ - co^ 



A direct use of Eq. (11), with |r| given as in Eq. (9), should yield back the 
phase 0(a;). As is well known one tacitly assumes that ln[T(Q)] is analytical 
on the upper complex Q-plane. This fact however is immediately granted by 
the fact that o is proportional to the imaginary part of the refractive index 
of the medium, and it is well known that the refractive index must obey 
Kramers-Kroning relation Eq. (11). We therefore used Eq. (11) in order to 
recover the phase of the transmittance. More specifically, we took advantage 
of a publicly available Matlab script [18] which serves exactly to that end. 
The final result in terms of modulus and phase of the transmissivity T is 
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Fig. 6. Gas cell dimensions. 



shown in Fig. 5. 



5 Physical analysis of the gas cell used in the monochromatization pro- 



cess 



An analysis of the gas cell used in the monochromatization process is par- 
ticularly important, since in the VUV range one cannot install a window 
separating the undulator pipe from the gas cell. The solution is based on 
the exploitation of a differential pumping scheme. The analysis presented 
here is based on [[l^-jlsl], which where used in the design of the LCLS gas 
attenuator. Fig. 6 shows the schematic of the gas cell that will be actually 
used in the monochromatization process. We will denote this cell as volume 
"0". 

Fig. 7 shows how the gas cell is integrated in the differential pumping 
system. Two differential pumping cells are locates symmetrically around 
voulme "0", after main cell (system symmetric in both sides). These cells 
will be denoted as volumes "1" and, similarly to volume "0", will have a 
length of 30 cm and a radius of 1 cm. They will be pumped by scroll vacuum 
pumps with a pumping speed of 30 1/s. These three cells (volumes "1" and 
"0") are installed by exploiting the transverse offset given by the magnetic 
chicane, which is limited to a maximum of 1 cm due to the requirement on 
the delay between electrons and X-rays. The system connected with vacuum 
chamber of magnetic chicane, which will be denoted as volume "2", through 
apertures of radius a = 0.5 mm. At both ends of the vacuum chamber, the 
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ISchematic ot the gas monochromator system 
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Fig. 7. Sketch of differential pumping system 

VUV are allowed to enter and exit through orifices with the same diameter. 
The dimensions of these apertures can obviously be changed. This means 
that if gas cell is not filled the orifices can be eliminated by fully opening 
them. When the gas cell operates, the orifices are in instead. The vacuum 
chamber volume "2" will be pumped by a turbo pump with a pumping 
speed of 300 1/s. As a result we have differential pumping system which 
includes two stages per each side: one of the two volume "1" cells as first 
stage, and the volume "2" vacuum chamber as second stage. A gas outflow 
of 10^^ atoms/s to the undulator vacuum chamber is obtainable, and is the 
same which was calculated for the LCLS gas attenuator, upon which the 
present analysis is based. 

A technical problem to be overcome pertains the diameter of the vacuum 
pump port. For different pressures different kind of pumps will be used. 
As said above, the two volume "1" cells will be at a pressure of 0.003 Torr. 
In this case scroll pumps will be used, with a port radius of 2 cm. Volume 
"2" will be at a pressure of 0.00001 Torr, in which case a turbo-pump will 
be employed, with a port radius of 5 cm. For the vacuum chamber in the 
magnetic chicane, volume "2", the turbo pump can be operated without 
issues, but for the scroll pumps in the two volume "1" cells, one has to deal 
with a small transverse offset in the centimeter size, while the port radius of 
the pump is about 2 cm. However, such constraint is only on the side where 
the beam offset is present. On the other side. Fig. 7, an adiabatic geometrical 
transformer can be used connecting the side of volume "1" (1 cm radius and 
30 cm long) to the 2 cm radius tube needed for the pump connection. 

Let us assume that the gas used in the monochromator is Helium. For Neon, 



15 



all issues associated with the gas monochromator become easier. In fact, due 
to the larger photo-absorption cross-section of Neon, the pressure in the cell 
can be an order of magnitude lower. Assuming that the gas in all systems 
is at roughly room temperature, the pressure P and the gas density Wq are 
related by 



no[cm-^] = 3.5 ■ 10^^P[Torr] . (12) 

Given the length / = 30 cm of volume "0", that is the absorption cell, one sees 
that a pressure Pq ~ 1 Torr is needed to reach the column density UqI = 10^^ 
cm"^ assumed in the previus sections. 

The collisional mean free path led is related to the density and to the elastic 
cross section Og by Zco/ = l/Cwo^e). Given the diameter of He atom da = 2.4A, 
the He-He elastic collision cross-section can be estimated as Oc ~ n{2daf ~ 
7 • 10"^^ cm^. It follows that, at the pressure of 1 Torr the mean free path is 
Icoi ~ 5 • 10"^ cm. This means that the gas flow is collisional, in the sense that 
the mean free path Icoi is much shorter than the orifice radius a. As a result, 
the gas exhaust through the orifices will occur in hydrodynamical manner 
(and not in Knudsen manner), and the number of atoms leaving the cell 
through each of the orifices can be evaluated as 



Qi ~ Q5na^nQSg , (13) 

where Sg is speed of sound, and the multiplicative factor 0.5 accounts to the 
reduction of the local speed of sound and density at the transmission point. 
Taking the density in volume "0" Wq ~ 3 • lO^^cw"^, ~ 3 • 10^ cm/s, and 
a = 0.5 mm one can calculate the throughput as Qi = 3 • 10^^ atoms/s. 

The gas density in volumes "1" depends on the pump speed of the corre- 
sponding pumps. The gas-load equation reads 



Qi = Spfii , (14) 

where Sp is the pumping speed in //s, and n\ is the gas density in volumes 
"1". Taking, as said above, the pumping speed equal to Sp = 30//s and the 
throughput Qi = 3 • 10^^ atoms/s, one can calculate the density and the 
pressure in cell "1" as as n\ ~ lO^^cm"^ and Pi ~ 3 • 10"^ Torr. 

Consider now the transition between cell "1" and volume "2". Here the 
mean-free path in the incoming gas is much larger than the orifice size 
a — 0.5 mm. The flow through the orifice under such circumstances is 
Knudsen flow. The vast majority of the atoms entering the volume "2" from 
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the cell "1" hits the walls of volume "2" and get scattered. These atoms are 
responsible for the density value n2. The gas outflow is 



Q2 ~ a'ni ^jlnkTIMue (15) 

Using parameters Ui ~ lO^^cm"'^, a = 0.5 mm, and T ~ 300 K, this yields 
Q2 ~ 5-10^^ atoms/s. We assume a pump speed for volume "2" of Sp ~ 3001/s, 
which can be easily reached with various types of pumps, if the surface area 
of the vacuum port in the side walls of volume is greater than 100 cm^. The 
density and pressure in volume "2" turns out to be then n2 ~ 2 • 10^^ cm"^, 
and P2 ~ 10"^ Torr respectively. 

We can now estimate the gas flow from the gas monochromator to the 
undulator beam line. The gas outflow to the undulator can be calculated in 
the same manner as the gas flow from the cell "1" to volume "2" and does 
not exceed 



Qun ~ a^n2 ^l2nkTIMHe ~ 10^^ atoms/s. (16) 

or equivalently Qun ~ 3 ■ 10"'' Torr 1/s. For a pumping speed in the undulator 
beamline Sp ~ 300Z/s, the final pressure will then be P^n ~ 10"^ Torr. 



6 Theoretical background of line profiles in the VUV absorption spectra 
of the rare gases 



In Section 3 we used the expression in Eq. (5) for our cross-section calcula- 
tions. In this Section we give a justification of that expression based on [24]. 
Let us consider for simplicity the Helium case below the second ionization 
threshold {Oa = 0). The PEL radiation can be treated as a time-dependent 
perturbation of a two-state system, the first state |1 > corresponding to the 
Helium ground state, and the second state |2 > corresponding to the final 
ionized state. In this case (see, any textbook of Quantum Mechanics, for 
example Izsjl ), we can make use of the Fermi's Golden Rule to calculate the 
transition rate between these states. Let us indicate with < 2|p|l > the matrix 
element of the perturbing potential, where the dipole moment operator p is 
given by 



V 



(17) 
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Tj being the transverse displacement of the jth electron in the atom. The 
transition rate can be found with the help of the Fermi's Golden Rule as 



W^2 = T^J-\<2\p\l>\'I{cv^2) 




(18) 



where I is the intensity, and cou is the transition frequency. A relation be- 
tween absorbed power P and incoming intensity I can be found by multi- 
plying both ends by the photon energy ha)u> yielding 



Since by definition the cross section is the ration between absorbed power 
P and incoming intensity I, it follows that the cross section is proportional 
to the square of the transmission amplitude. In order to calculate the cross- 
section in Eq. (5) we therefore need to calculate the square of the transmission 
amplitude. 

The interaction with the electromagnetic pulse has the effect of ionizing 
the atom. As discussed in Section 3, at some particular energies, ionization 
can happen either directly: He + hco — > He'^ + e~ or through a discrete 
autoionization state He' : He + hco — > He' — > He'^ -I- e". In the example case 
of Helium below the second ionization threshold considered here, such 
particular state He' corresponds to a double electronic excitation where the 
two electron principal numbers become ni = 2 and n2 = 2, 3.... 

In the quantum-mechanical language, what was just said means that the 
electromagnetic perturbation couples the ground state |i/^o > to the quasi- 
continuum states 3 {|^ >} of the unperturbed Hamiltonian Hq either directly, 
or through another discrete eigenvector of Hq, the autoionization state, 
which will be indicated with \(p >. The autoionization state \(p > is cou- 
pled to \k > through the coupling potential W, whose explicit determination 
will not concern us here and is in fact the subject for a density functional 
theory approach. We define the coupling strength between \(p > and \k > 
through Wk =< k\W\(p >. The coupling between \(p > and \k > through W 
modifies the structure of the quasi-conttnuum states \k >.ln fact, given the 
total Hamiltonian H = Hq + W, we can indicate its eigenvalues and eigen- 
vectors respectively with E^, and li/^^, >, and find the components of \ipi_i > 

^ Following ||2^ we discretize the continuum states of Hq by replacing them with 
discrete states ]k >, with eigenvalues Ej^, separated by a certain energy 6. When 
6 — > one obtains back physical results. 




(19) 



18 



on \(p > and \k > as fl^ 



< (pl^P^i >-- 



< k\^^, >-- 



1 + [w- 



-Ek' 

w/{E,, - Ek) 



1/2 



1 + 



-.1/2 ' 



(20) 



together with the eigenvalue equation 



E 



(21) 



In order to obtain Eq. (20) and Eq. (21), we assumed with [24(] that all Wk are 
equal and real, i.e. Wk = w, that < (/)|Ho|(/) >=< >=< k\]N\k' >= 0, and 

that < k\Ho\k >= Ek = kb. Equivalently, always following uM> we can write 
Eq. (20) as 



w 



< k\ip^ >-- 



+ (f )' + E 
w^/{E^, - Ek) 

w2 + (f )^ + E 



ll/2 



1/2 ' 



(22) 



where the transition rate F is F = 2nw^/{h6), and in the limit for 5 — > 0, 
up-fd is constant and equal to hT/{2n). 

For each choice of /i we are left with a two-level system with |1 >= \\po > 
and |2 >= >. The transition amplitude between state |1 > and state |2 > 
can be written as 



< 2|p|l >=< ip^,\p\ipo > 

=< ip^\(p >< (p\p\\po >+Yj^ '^''I^ ^IPl'^o ^ ' ^^^^ 

k 

which can be re-written with the help of Eq. (22) as 
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< 2\p\l >= 



+ [ff + E 



1/2 



(24) 



Finally, the use of the eigenvalue equation Eq. (21), the definition v -< 
(p\p\ipo > and the simplifying assumption that <k\p\ipo >= v' independently 
of |A: > allows to write 



< l\p\\ >= 



vw + v'E, 



+ (f )^ + £ 



1/2 



(25) 



Taking the limit for 5 — > 0, the quasi-continuum > becomes a contin- 
uum of states, and with the help of the reduced variables & = E^I{hYll) and 
c\ = 6v/{nzvv') = 2wv/{hTv') we obtain 



I < 2\p\l > P \q + £p 

= TTW ^'^^ 

In Eq. (26), E^, which is the energy of the continuum state, and can be written 
as the difference Epu - £12, with £12 = hco^ being the a theoretically defined 
resonance energy, and Eph being the photon energy of the electromagnetic 
pulse. By this, £^ can be positive or negative. This is consistent with previous 
choice to take the discrete energy level \(p > as the origin of the energy axis, 
that is < (p\Ho\(p >= 0. 

We reproduced the well-known result that the probability of exciting the 
atom from the ground state follows a Fano profile as a function of Epu. As 
said before, the cross-section is proportional to the square of the transition 
amplitude, and can therefore be written as 



o = (^b j-^ , (27) 

which corresponds to Eq. (5) with = 0. A generalization to the case with 
Oa non zero can be obtained considering that transitions to the continuum 
may not interact with the autoionization state, as discussed before. This can 
be readily done by introducing a the background cross-section Oa, and thus 
recovering Eq. (5). 
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7 Semi classical treatment of light propagation through a cell containing 
resonantly absorbing gas 



In this article we model the interaction of the gas in the cell with the FEL ra- 
diation response as the response of a macroscopic medium under the action 
of classical fields. Neglecting magnetic interaction and non-linear effects, the 
electric field E of the incident pulse induces a macroscopic polarization vec- 
tor P in the medium, which is linearly related to E through the susceptibility 
X according to the well known relation 



P = eoxE . (28) 

The susceptibility ;^ is a function of the the frequency of the applied field, 
and is related to the complex refractive index n via 



n^ = l+X. (29) 
and to the dielectric constant through 

e = eo(l + x) ■ (30) 

The determination of the susceptibility x depends on the dynamics of the 
system. For example, as is well-known, if the perturbation excites transitions 
from the ground state |1 > to an excited state |2 > of the atoms in the 
medium, one finds quantum-mechanically the usual result in agreement 
with the damped oscillator model: 



N\ < 2\p\l > p / 1 \ 

Xico) = — — 31) 

SeohV \coi2 - CO - ly ) 

where N is the number of atoms in the volume V, the operator p has been 
defined in Eq. (17), a»i2 is the transition frequency, and y models the radiative 
loss rate. Eq. (31) obeys, as well-known, a complex Lorentian relation. 

In principle one may calculate an analogous of Eq. (31) in the case when |2 > 
is a state in the continuum, coupled with atom in a doubly excited state, that 
is our autoionization state. In this case there are two contributions to the 
susceptibility. One arises from the dipole moment between the two states 
of the atom, and the other from the response to direct photo effect, that 
is the electromagnetic transition from ground state to continuum. In this 
article we do not present this direct calculation, and we make use of general 
properties of the susceptibility instead. 
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As is also well known, the susceptibility Xi and as a consequence the refrac- 
tive index and the dielectric constant, obey Kramers-Kroning relations, and 
it is always possible to recover the imaginary part of one of these quanti- 
ties from the a-priori knowledge of the real part, or viceversa. The typical 
example is the relation between real and imaginary part of Eq. (31). In our 
case, the a-priori information available is the knowledge of the cross-section, 
which obeys a Fano profile. The cross-section and the imaginary part of the 
refractive index are related through: 

Im[fcn] = ycT (32) 

where no is the gas density and = cu/c. As a result, as explained in Section 4, 
the prolongation of the logarithm of the transmittance to the upper complex 
plane is analytical, and the phase of the transmittance can be recovered 
from the knowledge of the cross-section. Once the transmittance T{oS) is 
known, the product with the incident field E(a;) yields the transmitted pulse 
Y.{oS)T{oS), which can then be characterized in the time domain by means of 
a Fourier transform. 



7.1 Collision line broadening 

As seen before, the knowledge of the cross-section allows for a full char- 
acterization of the medium behavior in the frequency domain in terms of 
transmittance. This is enough for our purposes, since the transmitted pulse 
can be recovered in the time domain by applying a Fourier transform. Simi- 
larly, the time-domain evolution of the macroscopic dipole moment P{t) can 
be found by Fourier transforming Eq. (28), which is based on the knowledge 
of x{(^)f recovered through the cross-section, and of the incident pulse in the 
frequency domain. The macroscopic dipole moment P{t) decays in a time 
scale of the order of ~ 100 fs, the inverse of the ionization rate. 

In this paragraph we will consider collision line broadening as a possible 
decay mechanism, and show that in our case of interest it can be neglected. 
Following, for example, ^2^, if the mean free path between two collision is 
much longer than the gas cell dimensions, one can consider the atoms in the 
gas moving along straight-line trajectories, colliding in a time which is much 
shorter than the radiative lifetime of the atomic excited state. Moreover, the 
time between different collision is long compared with such lifetime too. 
Then, the collision can be modeled as a loss of coherence due to a sudden 
phase jump of the atomic excited wavefunction. The time Tcoi between two 
collisions, i.e. the time between sudden phase jumps, is the inverse of the 
coUisional rate y'^^ Always following [26] we demonstrate that collision line 
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self -seeding scheme with gas monochromator for the LCLS-II basehne soft X-ray undulator 



electron beam 
energy 3.5 GeV 
300 fs 
0.25 nC 



tunable gap beta function 10 m 

undulator cell length 3,7 m 

period 60 mm undulator segment 

K = 1.3-12 length 3 m 



phase I 
22 cells 



SASE 




undulator 


3.7 m 


7 cells 




19,2 nm 


gas cell 


K-7.9 


with He 




1 torr 



6 cells 9 cells 



I 



19.2 nm 4.8 nm 
K = 3.8 



self -seeded soft 
X-ray pulse 
at 4.8 nm 
energy 1 mJ 
bandwidth 0.01 1 



Fig. 8. Design of the LCLS-II baseline imdulator system for generating of generating 
of highly monochromatic, high power soft X-ray pulses. The method exploits a 
combination of a self-seeding scheme with gas monochromator and gap-timable 
undulator for harmonic generation. 

broadening can be neglected in our case. In order to do so, we first observe 
that, for hard-sphere collisions between atoms of mass M and radius and 
assuming a gas density Wq/ the inverse of the coUisional rate is 



7col = col ^ 



SRjtio 



(33) 



In our case for the Helium gas we have Ra = 1.2 A, Uq = 3 ■ lQ^^cm~^ (for a 
pressure of 1 Torr at room temperature) Mc^ — 410^ eV, kT — 0.03 eV, and as 
a result we obtain Xcoi - 400 ns. Since this time scale is much longer than any 
other time scale in our system, we can neglect this effect. In other words, the 
decay of the macroscopic polarization is due to spontaneous radiation only. 



Operation of the LCLS-II soft X-ray FEL by use of the self-seeding 
scheme with gas monochromator. A feasibility study. 



The autoionizing resonance that we consider here as an example is in the 
20 nm-range. In order to generate shorter wavelengths we need to go to 
higher harmonic numbers. A possible way of doing so is to use a two-stage 
output undulator, with the second stage resonant to one of the harmonics of 
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magnetic chicane dimensions for 
gas monochromator at LCLS-II 




Fig. 9. Magnetic chicane dimensions for gas monochromator installation at LCLS-II 

the first one. Fig. 8. The magnetic chicane dimensions are sketched in Fig. 9. 
In our feasibility study we will consider, for exemplification purposes, the 
LCLS-II setup equipped with a 80 m long gap-tunable undulator, which can 
be operated within the photon energy range between 30 eV to 1100 eV at 
3.5 GeV electron beam energy. Due to specific of Front End Enclosure (FEE) 
system of LCLS-II will be possible direct the radiation into experimental 
stations only starting, around the 200 eV photon energy range. With the 
chosen parameters [£J (see Table 3 and Fig. 10), the amplitude of the 4th 
harmonic of the density modulation at the exit of first stage of the output 
undulator is high enough and dominates significantly over the amplitude 
of shot noise harmonic. This modulation density serves as an input signal 
for the second part of output undulator, which is indeed resonant with 
the 4th harmonic. An important feature of our design is that no dispersion 
section is introduced between the two stages in the output undulator. The 
advantage of this method of harmonic generation in self-seeded FELs is in 
the simple hardware required (only a short magnetic chicane for the gas cell 
installation). However, a gap-tunable baseline undulator is needed. 

Here we focus on the implementation of our self-seeding scheme exploiting 
a wake monochromator to generate fully coherent soft X-rays for LCLS- 
II. The tentative design for the gas monochromator self-seeding technique 
at the LCLS-II aims at generating fully coherent radiation at about 5 nm. 
The main parameters for the nominal 5 nm case are listed in Table 3 in 
the next Section. Full 3D have been performed in order to confirm the 

^ We used a not updated set for the LCLS-II SXR undulator parameters for our 
simulations, but the same wavelength is achievable for the current design too. 
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Table 3 

Parameters for the nominal pulse mode of operation used in this paper. 

Units 



Undulator period 


mm 


60 


K parameter 




1.3-12 


Wavelength (fundamental) 


nm 


19.2 


Energy 


GeV 


3.5 


Charge 


nC 


0.25 







Fig. 10. Electron beam characteristics at the entrance of the setup at the LCLS. (upper 
left) Current profile, (upper right) horizontal and vertical geometrical emittance. 
(lower left) energy profile, (lower right) rms energy spread profile. 

scheme feasibility. The simulations is performed with the code GENESIS 
1.3 code jz^I, which uses as input the beam parameters obtained in start to 
end simulations | l28l ]. The beam distributions are shown in Fig. 10 showing 
the current profile, the horizontal and vertical geometrical emittances, the 
energy profile and the rms energy spread profile. The average betatron 
function used was jS = 10 m. 

In all simulations we assumed that the gas used in the monochromator is 
Helium, and that the energy of the VUV EEL photons is 64.4 eV, or 19.2 
nm, corresponding to the n = 4 autoionization profiles for the {sp,2n+) 
autoionizing series of Helium, discussed in Section 3. 
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9,0 19,1 19,2 19.3 19.4 19.5 




X[nm] 



s[pm] 



Fig. 11. Pulse spectrum (left plot) and power (right plot) before the gas cell. Grey 
lines refer to single shot realizations, the black line refers to an average over one 
hundred realizations. 




X[nm] 



s\pm] 



Fig. 12. Pulse spectrum after the gas cell (left plot) and seed power (right plot). 
Grey lines refer to single shot realizations, the black line refers to an average over 
one hundred realizations. 

The spatial and spectral structure of the LCLS-II pulse at 19.2 nm wavelength 
and at the location of the gas monochromator are showr{£] in Fig. 11. The 
effect of the passage of the photon pulse through the gas cell is shown in Fig. 
12, where one can clearly see the bandstop filtering in terms of spectrum, 
and its counterpart in the time domain in terms od seeding pulse. The FEL 
radiation pulse has a transverse FWHM size of about 0.5 mm at the exit of 
the differential pumping system, the red line in Fig. 13. 

While the photon beam passes through the gas cell, the electron crosses the 
chicane, where the microbunching is removed. In fact, for Gaussian local 



^ From Fig. 11 and Fig. 13 we can estimate a photon density of 10^^ cm~-^. It 
should be noted that the resonance line is about two orders of magnitude narrower 
(10~^) compared with the SASE bandwidth (5 • 10~^), while the maximum cross 
section reaches about 10 Mb. As a result, attenuation takes place mainly due to the 
background cross-section, of about 1.5 Mb. Recalling that 1 Mb is 10~^^ cm^ one 
can confirm the validity of the perturbation theory approach used to calculate the 
overall transmission function of the gas cell. 
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Fig. 13. Output size of the FEL radiation pulse. The black line is the photon beam 
profile before the monochromator, the red line refers to the photon beam profile 
after 4 m propagation, that is at the exit of the differential pumping system. 

energy spread, the amplitude of the density modulation a at the chicane exit 
is given by 



where Uq is the amplitude of the density modulation at the entrance of 
the chicane, R56 is the momentum compaction factor and A is the reduced 
wavelength. In our case, parameters of interest are the dispersion R55 = 
Lzy0^ ~ 0.2 mm (insuring a delay equal to 0.1 mm), and the relative energy 
spread Ay/y ~ 0.03%, corresponding to an energy spread of about 1.5 MeV 
at an energy of 3.5 GeV. For a wavelength of 19 nm, these parameters lead 
to the suppression of the beam modulation by a factor of about exp(-200). 

Following the chicane, the electron beam is recombined with the photon 
beam at the entrance of the output undulator. The first stage of the output 
undulator, consisting of the first 6 cells, is tuned at the fundamental, while 
the second stage, consisting of 9 cells, is tuned at the fourth harmonic. Fig. 
8. The output characteristics in terms of power and spectra at the exit of the 
first stage are shown in Fig. 14. 



a = Uq exp - 



2 72 ;t2 



(34) 
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x[nm] s[pm] 

Fig. 14. Output characteristics for the second imdulator tuned at the fundamental 
harmonic. Spectrum (left plot) and power (right plot). Grey lines refer to single shot 
realizations, the black line refers to an average over one hundred realizations. 




Mnm] s[ijm] 

Fig. 15. Output characteristics for the second imdulator tuned at the fourth har- 
monic. Spectrum (left plot) and power (right plot). Grey lines refer to single shot 
realizations, the black line refers to an average over one himdred realizations. 

The final output at the fourth harmonic is shown in Fig. 15, showing power 
and spectrum, and in Fig. 16, showing energy and variance of the electron 
energy fluctuations. A mj, fully coherent pulse with a bandwidth in the 
order of the 0.01% is produced at a wavelength of 4.8 nm. 



9 Conclusions 

In this paper we propose a method for controlling the line width of VUV 
and soft X-ray SASE FELs that offers simplicity and flexibility, and can be 
added to the baseline undulators of many facilities without significant cost 
or design changes. Monochromatization down to the Fourier transform limit 
of the radiation pulse can be performed by exploiting an almost trivial setup 
composed of as few as two components. The key components of such scheme 
include only a cell containing dilute noble gas, and short magnetic chicane. A 
great advantage of our method is that in includes no path delay of radiation 
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1.2x10^- 




Fig. 16. Output characteristics for the second undulator tuned at the fourth har- 
monic, as a function of the length of the second undulator. Energy (left plot) and 
variance of the energy fluctuations (right plot). Grey lines on the left plot refer 
to single shot realizations, the black line refers to an average over one hundred 
realizations. 

pulse in the monochromator. Implementation of the proposed technique 
in baseline undulator will not perturb the baseline mode of operation. We 
present an illustration of the scheme for the LCLS-II soft X-ray beam line, 
although other facilities like FLASH-II, SwissFEL, FERMI and SPARC may 
also benefit from it. 
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